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Quantum decoherence in a hybrid atom-optical system of a one-dimensional 
coupled-resonator waveguide and an atom 
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Decoherence for a one-dimensional coupled-resonator waveguide with a two-level system inside 
one of resonators, induced by their interaction with corresponding environments, is investigated. 
Each environment is modeled as a continuum of harmonic oscillators. By finding the eigenstates of 
the hybrid system, which is the dressed state of the hybrid system, we calculate the lifetime of one 
excitation, which characterizes the existence of quantum coherence in such hybrid system and the 
basic quantum nature. 
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I. INTRODUCTION 

Any discrete state coupled to continua of states is sub- 
ject to decay. The intrinsic dynamics of such quantum 
system is irreversible, i.e. once a system is initially in the 
discrete state, it never returns to the initial state spon- 
taneously. Such kind of phenomenon is also described 
by the so-called resonant tunneling process. There are 
many examples of such processes in different branches of 
physics. A typical example is an excited two-level system 
interacting with the modes of electromagnetic fields. 

Spontaneous emission is normally regarded as a loss 
and a decoherence mechanism. Recent works show that 
spontaneous emission of the excited two-level system can 
be exploited to influence the coherent transport proper- 
ties of single photon in a one-dimensional (ID) waveguide 
due to the interference between the spontaneous emission 
from two-level systems and the propagating modes in the 
ID continuum [l[. For a discrete system, the simplest 
model which possesses the resonant tunneling process is 
the so-called Anderson- Fano-Lee model [2|-l5|, in which 
the continuum is formed by a linear chain of sites with 
the nearest-neighbor interaction and the rate of emission 
is modified due to the change of the density of state. 

With the development of microfabrication technology, 
the platforms, such as defect resonators in photonic crys- 
tals |7H10| and coupled superconducting transmission line 
resonators [ll| - [l4| , are promising candidates for realizing 
a waveguide (or an array) [g], under the tight-binding 
approximation. Since a ID coupled-resonator waveguide 
(CRW) possesses a band-gap spectrum and can trans- 
mit a wave packet of hght, a single-photon quantum 
switch [1^ [lg| , made of a controllable two- level or three- 
level system, has been studied using a discrete-coordinate 
approach. It is then found that the transmission of a 
single-photon in a ID CRW can be switched on and off 
by modulating the energy-level spacing of the two-level 



'Corresponding author; Electronic address; |zzhoulan@gmail.com| 



system (TLS) with a high-frequency signal |17| . 

The system with TLSs inside a ID CRW [IMl, is 
investigated under the condition of ideal resonators and 
ideal couplings of the resonator modes to the respective 
atomic transitions. Such a system is closed, and the exis- 
tence of superpositions prescribed by quantum mechanics 
is valid. All the results in Refs. [l^ [la] are obtained by 
seeking the stationary states of the system with a TLS 
inside ID CRW. However a realistic quantum system can 
rarely be isolated from its surrounding completely, rather 
it is usually coupled to the external environment (also 
called "heat bath" ) with a large number of degrees of 
freedom. There are two main loss processes which are 
serious obstacles against the preservation of quantum co- 
herence over long period of time: spontaneous emission 
from the excited state to the ground state due to its inter- 
action with the modes outside the resonator, and leaking 
out of photons of the resonator mode. Obviously, when 
the external environment is taken into account, the deco- 
herence of every resonator and the TLS would result in 
the incoherent or dissipative propagation of the incident 
photon. 

It is pointed in Refs. |15l-[l7{ that the decoherence or 
dissipation can be divided into two categories: one in- 
fiuences the free propagation of the single photon, and 
another infiuences the scattering process which broadens 
the line width. 

Therefore in the present work, we shall investigate the 
infiuence of environment to the decoherence of ID cou- 
pled resonator with a two-level system inside. We mainly 
focus on two issues, one is the lifetime of each eigenstates 
characterizing the time-scale of quantum coherence in a 
ID CRW with a TLS inside, and another is the the rate 
of decay through introducing the leakage rate in each res- 
onator and the decay rate of the TLS, which influences 
the free propagation of the single photon. The reflection 
spectrum is also obtained and it is found that the dis- 
sipation lowers the peak of the resonance and broadens 
the line width. However, the total reflection can still be 
achieved when the leakage rate in each resonator is equal 
to the decay rate of the TLS. 



This paper is organized as follows. In Sec. HIl we intro- 
duce a microscopic model where both the resonators and 
the TLS are coupled to its own surrounding through an 
exchange interaction. In Sec. IIII[ we briefly review the 
stationary states of a CRW with a TLS inside in Ref. [l^ . 
In Sec. IIV) we derive the characteristic time for a single- 
photon staying in the system with a TLS embedded in 
ID CRW. In Sec. |Vl we investigate the impact of dissi- 
pation on scattering amplitude. Conclusions are made in 
Sec. ED 
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FigureH] shows a schematic diagram of what we con- 
sider in this paper. The total system includes the system 
S and the environment B. The environment B doesn't ex- 
plicitly show in FigUl The system S refers to a ID CRW 
with a TLS inside one of the resonators as discussed in 
Ref.[la| and the environment B refers to all the subsys- 
tems interacting with each resonator in the ID CRW and 
the TLS. In FiglHthe wavy line with an arrow indicates 
that a bath interacts with this subsystem. 

Each resonator in the CRW is modeled as a harmonic 
oscillator mode with frequency uJc- Due to the overlap 
of the spatial profile of the resonator modes, photons 
can hop between neighboring resonators. Introducing the 
creation and annihilation operators aj and aj for jth res- 
onator, the Hamiltonian for the CRW is given by 
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where the intercavity coupling constant ^ is the same 
for all neighboring cavity-cavity interactions. A TLS 
with transition frequency ft is located in the 0th res- 
onator. The atom-resonator interaction is described by 
the Jaynes-Cummings model 



i7i-r!|e)(e| + J(|e)(.g|ao+H.c.) 



(2) 



under the dipole and rotating wave approximations. Here 
|e) and l^) are the excited state and ground state of the 
TLS, respectively. The dynamics of the system S is gov- 
erned by the Hamiltonian 



^s — He + Hi. 
By employing the Fourier transformation 
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FIG. 1: (Color online) Schematic illustration of the model, 
where a two level system (a) is inside a ID coupled resonator 
waveguide (b). The two level system and each of resonators 
are coupled to its corresponding environment, which is mod- 
eled as an infinite set of harmonic oscillators. 



where the lattice constant is assumed to be unity and 
fife = ojc — 2^ cos k is the well-known Bloch dispersion 
relation. Here, the Hamiltonian He given by a diagonal 
matrix describes the extended states of the continuum. 
The third term in Eq. ([S]) is responsible for the interaction 
between the TLS and the continuum. The Hamiltonian 
Hs describes that a quasiexcitation is created or anni- 
hilated in the fc-th mode, and the TLS transits from its 
excited state to the ground state or vice versa. Obviously, 
it is easy for the TLS to transit to its ground state, and 
very difficult to go back to its excited state. Therefore, 
the CRW with a TLS inside is a typical system of quan- 
tum dissipation phenomenon. 

The environment B is modeled as a set of infinite num- 
ber of harmonic oscillators. In this paper, we assume 
each resonator of the CRW is coupled to an individual 
bath of harmonic oscillators. The dynamic of the jth 
resonator and its corresponding environment is governed 
by the Hamiltonian 
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q[j] 



H.c. 



(6) 



where bqy] and b ,-, are the annihilation and creation 
operators for the gth bath oscillator attached to the jth 
resonator of the CRW, ujqy^ is its frequency, and gg[j] the 
coupling strength. The total Hamiltonian describing the 
interaction between the CRW and environments reads 



Hec - X! ^ER- 



(7) 



Hamiltonian Hs is transformed into a k-space represen- 
tation 

i?s = E ^kaiak + n |e) {e\ + ^=Y^ (|e) {g\ Uk + H.c.) 
k ^ k 

(5) 



We also introduce the exchange interaction between the 
TLS and its environment, which is described by the 
Hamiltonian 



Hea = ^ i'qd\dq 



^/3„(4|5)(e| + |e)(5|d„), (8) 



where d„ and dj^ are the annihilation and creation oper- 
ators corresponding to the qih mode with frequency i^q. 
The Hamiltonian of the total system S+B reads 



^SB — Hs + H] 



EC 



Hea- 



(9) 



Although there exists energy exchange between the sys- 
tem S and the environment B, the total number of exci- 
tations is preserved. 



III. STATIONARY STATES OF A CRW WITH A 
TLS INSIDE 

In the one-excitation subspace, the stationary states 
of system S are either the localized states around the 
location of the TLS or a superposition of extended prop- 
agating Bloch waves incident reflected and transmitted 
by the TLS embedded in the CRW. Assume a photon is 
coming from the left of the TLS with energy 



ni 



ujc — 2^cosfc. 



The stationary state of the system is then 



m, 
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Uk{j)a] \0g) +Uek |0e) , 



(10) 



(11) 



where |0) is the vacuum state of the cavity field, u^k is the 
probability amplitude with the TLS in the excited state 
and no photon, and Uk (j) is the probability amplitude 
for the TLS to move to the lower state, emitting a photon 
into a mode of the jth resonator. From the eigenvalue 
equation Hg \i^k) = ^k \^k), together with the bosonic 
commutation relations, we derive a system of coupled 
equations among the amplitudes 

(Vk + JGk6,o) Uk (j) = e K(J + l)+uk (j - lE2a) 

Uek = GfeWfe(O), (12b) 

with the Green function Gk = Gk{^) — J / {^k — ^) and 
Vk = ujc ~ ^k- Equation (|lll) has the solution in terms 
of incoming and outgoing waves with amplitudes 



Uk (j) 
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The reflection and transmission amplitudes yield 



rk 



Sk 



Sk - 1, 

2i^ (fifc — Q) sinfc 
2i^ (tlfc - f}) sin fc - J2 ■ 



(13) 

(14a) 
(14b) 



Eq. ([T2|) implies that the transmission coefficient is sen- 
sitive to transition frequency of the TLS, as the forward 
and backward propagating modes within the CRW are 
coupled via the TLS. 

The periodicity of the CRW in r-space gives arise to 
a continuum with Bloch waves grouped in energy bands, 
which is broken due to the coupling between the TLS 



and the CRW. Local modes are produced, and their cor- 
responding eigenenergy is outside the energy bands. We 
denote the eigenvalue of the bound state as fi^, i.e. re- 
place the index k with k, so does the equation (jlip . 
Noticing that the probability of of bound state at infinity 
in coordinate space is zero, i.e. bound state is spatially 
localized, we can assume that the bound states of even 
parity have the following amplitudes 



"k (j) 






where the normalization constant C is 



C = 



tanh k 



J2 



(Ok - ^Y 



-1/2 



(15) 



(16) 



The bound state energy lies either below or above the 
continuum with the magnitude 



n^ =wc-2^e"''coshK, 



(17) 



where the value of k is determined by the following con- 
dition 

j2 ^ 2^6™'' {^-^k) sinh K. (18) 

The bound states exist when Ea.([T^ has solution. 



IV. SINGLE-PARTICLE 
SPONTANEOUS-EMISSION 

Quantum states inevitably decay with time into a 
probabilistic mixture of classical states due to their in- 
teraction with the environment, which comprises much 
larger systems or ensemble of states. In this section we 
derive the relaxation time for single combined photonic- 
atomic excitation in the system S by incorporating the 
effects of interaction with the environment B. Notice that 
the number of excitations is conserved in the total system 
(S-l-B). Therefore, in one-excitation subspace, two situa- 
tions will occur in the system S, one is that the single- 
excitation is in state \Q,k) or \^k)] and another is that 
no excitation is found in the system S, i.e. the TLS is 
in the ground state and no photon is in the resonator. 
For the first case, there is no excitation in the environ- 
ment B, and we use the notation jilfcOO) to describe the 
state of the total system, where the first in the ket 
indicates that the radiation modes coupled to the CRW 
are in the vacuum state, and the second in the |rifcOO) 
denotes the vacuum state of the harmonic oscillators at- 
tached to the TLS. When the system S is in \{)g) state, 
the environment B contains the single-excitation. Here, 
djj |G) (|G) = lOgOO)) represents that the single-particle 
has moved to the nth bath oscillator coupled to the TLS; 
and 6 , , \G) describes that the excitation in S is in the 
gth mode of environment B attached to the jth resonator 
of the CRW. Obviously, states ||f^feOO) , d^ \G) , b\^,^ \G)\ 



provide a complete basis and thus we can expand the 
wave function of the S+B at arbitrary time t in terms of 
this basis as 

l*W> = Y.^kit)\nkOO) + J2c„{t)di\G) (19) 

k n 

where Uk, Bgy^, and C„ are the time-dependent proba- 
bihty amphtudes for finding the entire system in its cor- 
responding states. 

Inserting the wavefunction (fT9)) into the Schrodinger 
equation with the governing Hamihonian in Eq. (|9]), we 
find a system of coupled linear differential equations for 
the amplitudes 

iilk = ^kUk + Yl 9qU]ul (j) B,[j] + Yl ulkf3„(Z]a) 



the exact integro-differential equation 



fc 

iCn = VnCn+'^PnUekUk, 



(20b) 



(20c) 



where the overdot indicates the derivative with respect 
to time. The evolution of Uk is coupled to S|j[j] and 
C„ via the coupling constant gq[j\u*j. (j) and u*i^$n- To 
remove the high-frequency effect, we make the following 
substitution 



Uk{t) = 4>k{t)e 
BqU]{t) = b'gy^{t)e 



— iflkt 



(21) 
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Then in the interaction picture, equation (|20p is rewritten 



as 



i<Pk = J2 9qU]^k U) b'q[jf '^^"^ +J2Kkl3nCn( 
jq n 
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where A|v = Wgui — Vtk, 5nk — Vn~^k are the detunings 
between the system S and the environment B. Formally 
integrating the equations for bqu^ and c„ in the above 
equations yields 

%m = -*5g[,]E^^(j) /'^'^■(^)^'^''^'^^ (22a) 
k •' 

Cn = -iY, PnUek I ^k (t) e'^^"^ dr. (22b) 

fc '' 

Inserting Egs. (I22p into the equation for 0^, we obtain 



kjq 
kq 



iA^ t 



/ </>fc 
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(r) e' '''''^dr 



</)fc(T)e"''-^dT, (23) 



where the coupling constants ggtji and /3g are assumed to 
be real. 

Suppose that at time t = there is no interaction 
between the system S and the bath B, and the S-l-B is in 
the state jri^OO), corresponding to a single excitation in 
the system S. Mathematically, solving Eq. ([23]) is to solve 
a initial value problem. It is well known that such an 
initial value problem can be solved by Laplace transform. 
Denoting the Laplace transform of </>„ (t) by </>„ (s) and 
taking into account of the initial conditions (j)k (0) = dnk , 
we have 



1 



where 
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(24) 



The roots of the denominator in Eq. (|24| can be split 
into a sum of the singular and principal value parts. The 
principal part is merely included into redefinition of the 
energy. The singular part gives the decay rate due to 
the coupling to the bath B. Under the Wigner-Weisskopf 
approximation |l9|, the system decay is dominantly ex- 
ponential with rate 

r„ = TT ^ |U„ (i)l' Aj (r!„) + IT \Uen? ^A [^n) ■ (25) 

This rate is proportional to the modulus square of the 
coupling between S and B. 
The functions in Ea.([25| 



A,(^) = E-4i'^' 



'^qm) 



q 

are called the reservoir response (memory) functions, 
which are the spectral densities of the states b ,-, \G) 



q[j] 



and dj \G) weighted with the coupling strengths g; 
and Pg respectively. The memory function A (w) char- 
acterizes the spectral shape of the reservoir. Since the 
states in the reservoir are very dense (continuum), one 
can replace the sums over q by integrals, for instance, 
J2q -^ 2V/{2T:f f q^ smed(j)dedq. Then the memory 
functions read 

Aj (w) -^ gfj] (uj) py] (uj) , 
Aa(w) ^ ;3^(w)pa(w), 
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FIG. 2: (Color online) The decay rate via wavenumber k 
starting from one eigenstate of the system S, with parameters 
g = 0.1, /3 = 0.4, J = 1.5 , (a) wc = 5, fi = 6, (b) cue =^^5, 
which are in units of ^. 



where g?-, (w) are the couphng constants to states 
6^r , |G), whose density of state is given by pyj (uj). Obvi- 
ously, the decay rate is determined by two factors, 1) the 
form of A (w); 2) the width of the reservoir, which depicts 
the overlap of the eigenfrequencies of the system and the 
spectral shape of the reservoir. Usually, it is assumed 
that the reservoir is spectrally flat and the frequencies 
of the system of interest is deeply embedded in the con- 
tinuum, therefore, the argument of the spectral density 
is replaced by the eigenfrequency of the system r2„ in 
Eq. (1251) . For the sake of simplicity, we further assume 
that the coupling constant ggij] is independent of the en- 
ergy ujq[j] as well as the index j, i.e. let A[j]((u;) = 5(0;), 
and let Aa(j-o) = (3a(j-o). The constant coupling gives us 
the explicit relationship between the decay rate and the 
wavenumber of the system S 



irg 



n{/3 



(26) 



where the normalization of the eigenstates |ri„) has been 
used in deriving above equation. Obviously, when /3a = 
g, the decay rate is a constant for any eigenstates of the 
system S. For a given wavenumber n, when Pa is smaller 
than g, the rate decreases, which opens up a possibil- 
ity of prolonging the lifetime of the single-excitation in 
the system S. In Fig[21 we plot the decay rate F^ as 
a function of the wavenumber k. It indicates that the 
widths of the eigestates in Eg. p^ are different, which is 
a periodical function of the wavenumber fc, and the de- 
tuning between the TLS and the resonator biases the 
symmetry of the lineshape. By substituting the con- 
dition in Eq. (fT8| into the expression of the probability 
amplitude Uek, the decay rate of a bound state has the 



1- (J^ + 2^2 sinh 2k) 



form F«, = 7rg2 + TT (^^ - g 

It indicates that the decay rate is an increasing func- 
tion of K, and when the resonator resonates with the 
TLS, r„ = TT (2g2 + ^^) /3. Since the value of the imag- 
inary wavenumber is determined for given parameters 
{J,£,,uJc' ~ ^}, the decay rate of the bound state is a 
given constant. 



V. DISSIPATION ON SCATTERING 
AMPLITUDE 

To see the impact of dissipation of the system on the 
transmission/reflection coefficients, we begin with Hamil- 
tonian ^ in the configuration space. Here, it is assumed 
that each bath attached to its resonator is identical and 
the coupling strength between the resonator and its cor- 
responding bath is independent of the location index j. 
The state at arbitrary time is a superposition of four 
parts: the photon at the jth cavity with atom in the 
ground state a, \G), no photon in all cavities with atom 
in the excited state (T+ |G), the photon in the gth mode 
of the bath attached to the jth resonator o',., \G) and 
that of the bath coupled to the TLS cj \G) 

q 

The Schrodinger equation results in a system of coupled 
linear differential equations for the amplitudes 



i<^jk = UJc^jk - ^ {^3 + lk + *i-lfe) + J'^ekSoj + 2^ Qq^^) 
i^ek = fi^efe + J$Ofc + E /^9^9 

ibq = l3q<^ek + VqDq, 

where the dot denotes a derivative with respect to time. 
Applying the Fourier transform and expressing Aq^j^ , Dq 
with <I>j7£ , $efc respectively, the equations of motion in the 
frequency domain are obtained in the reduced dimension- 
ality 



{E -uc ~ SE) xpjk = -S, {ipj+ik + -ipj-ik)+- 



E-n- AE' 
(29) 
where E is the eigenenergy of the single-photon wave in 
the whole system. The eigenfrequency of the jth res- 
onator of the ideal CRW is renormalized into cue + SE, 
and the transition energy of the TLS is also renormalized 
into n + AE 



se=j: 



E-uj, 



., AE = y^ 



E-i^„ 



(30) 



which is the influence of the bathes on the state of each 
resonator and the TLS. The singular part of 6E and AE 
yields the dissipation factors 7c — irg^ (loc) Pc i^c) and 
J A — T^P (^) Pa {^)- The real part of 5E and AE con- 
tributes to the Lamb shift of the levels and change of 51, 
respectively, which is merely included in the redefinition 
of the energy. Then Eq. (P^ is reduced to the following 



set of equations 
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'j-lk, 



j'^tpOkSoj 



(31) 
We begin our analysis from the case with coupHng 
strength J — 0. In this case, Eq.((29|) has a complex 
dispersion curve that can be identified with the plane- 
wave solution $jfc (X e*('^-'~^*\ where E = E^ + iEi is the 
complex energy and k is the wave momentum inside the 
CRW. The real and imaginary parts of energy E satisfy 
Ej. = J7fc, Ei = 7c- When only one resonator is initially 
excited |^ (0)) — \n) (say n = 0), where \n) is the Wan- 
nier state localized at the site n, the field profile at time 
t is given by 



l*W) 



—iujct—lct 



^J5(2^i)e''?'^/2g-»/cg|^^ (32) 



kq 



where \k) is the Bloch state. The amplitude at the site I 
reads 



$ 



Ik 



-''''* Ji{2^t)e'^'^^^ 



(33) 



where Ji (x) is a Bessel function of the first kind of in- 
teger order I. Obviously, leakage rate in each resonator 
7c influences the free propagation of the single photon. 
The distance that the photon travels along the ID res- 
onator waveguide is depicted by the product of the group 
velocity and 7~^. 

We now consider the case with coupling strength J ^ 
0. Due to the coupling between the atom and the 0th 
resonator, a complex (5-like potential stands in the way 
that single photon travels. Consequently, the photon ex- 
periences scattering. We assume that the single-photon 
has momentum k initially. Within the allowed distance 
of photon traveling along the ID CRW, the transmission 
and reflection amplitudes, t and r, can be deflned via the 
asymptotes of the wave function 



tik = 



;«fci +j.fce-*'=^', j <0 



Sue 



ikj 



j>0. 



(34) 



By substituting the asymptotes of the wave function into 
Equation (|5T|) for j = ±1, sites, the reflection amplitude 
is obtained as 



.P 



rk 



2i^ sin k {ilk — fl) — {jA — 7c) 2^ sin k ~ J-^ 



(35) 



Equation (|35|) shows that as long as the rate of decay to 
each resonator 7c is equal to the decay rate of the TLS 
7^, a resonant scattering occurs when the incident energy 
of the single-photon is equal to 51, the transition energy 
of the TLS, i.e. when the single-photon incident from 
the left encounters the TLS, it is completely reflected 
back to the left, as shown via the green dotted line in 
Fig. 131 For unequal 7c and 7a, the decay rates lower the 
peak of the resonance, and the width of the line shape 
is broadened from J^/(2^sinfc) for ideal system S to 
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FIG. 3: (Color online) The reflection coeflicient versus the 
wavenumber k . The parameters are in units of ^ and are set 
as follows: 7c = 0.1, 7a = 0.4, J — 0.8, ujc — i^ — 5 for blue 
solid line; 7c — 0.1, 7a = 0.4, J = 1.5, ujc = SI = 5 for red 
dot-dashed line; 7c = 0.1, 7a = 0.4, J = 1.5, Uc = 5, SI = 6 for 
black dashed line; 7c = 0.1, 7a ~ 0.1, J = 1.5, Uc = 5, S7 = 6 
for green dotted line. 



[ilA - 7c) 2C sin k + J"^]/ (2^ sin k). In FigH we plot the 
reflection coefficient as a function of the wavenumber k. 
Comparing the blue solid line and the red dash-dot line 
with the black dashed line and the green dotted line, 
it can be found that the symmetry of the lineshape is 
determined by the transition energy SI of the TLS. As the 
coupling strength J between the 0th resonator and the 
TLS goes to infinity, the reflection coefficient approaches 
to one. 



VI. CONCLUSION 

In summary, we have studied the dissipative process of 
the CRW with a TLS inside, originating from the cou- 
pling to the environment. Our discussions are based on 
a simple decoherence scenario in which each resonator 
and the TLS individually interact with its own environ- 
ment. Each environment is modeled as a continuum of 
harmonic oscillators and assumed to be on the vacuum 
state initially. Since the coupling with the environment 
is generally weak compared to the system of interest, the 
Wigner-Weisskopf approximation applies, and the life- 
time of an excitation in the system S is thus obtained, 
which represents the time-scale of the transition from 
quantum to classical behavior. We further investigate 
the impact of the dissipation on the transport property 
of single-photon along the ID CRW for identical bathes 
attached to the CRW. The dissipation lowers the peak of 
the resonance and broadens the width of the line shape of 
the reflection spectrum except the case that magnitude 
of the leakage rate is equal to the decay rate of the TLS. 
The leakage rate in each resonator not only influences the 
free propagation of the single photon, but also results in 
the inelastic scattering of the single-photon together with 
the dissipation of the TLS. 
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